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THIRD QUARTER SECOND YEAR PROGRESS REPORT 

This report covers technical progress during the third quarter of the second year of NASA 
Sun-Earth Connections Theory Program (SECTP) contract “The Structure and Dynamics of 
the Solar Corona and Inner Heliosphere,” NAS5-99188, between NASA and Science 
Applications International Corporation, and covers the period February 16, 2001 to May 15, 
2001 . Under this contract SAIC and the University of California, Irvine (UCI) have conducted 
research into theoretical modeling of active regions, the solar corona, and the inner heliosphere, 
using the MHD model. 


SUMMARY OF ACCOMPLISHMENTS 

In this report we summarize the accomplishments made by our group during the first 
seven quarters of our Sun-Earth Connection Theory Program contract. The descriptions are 
intended to illustrate our principal results. A full account can be found in the referenced 
publications. 

1.1. Modeling the Large-Scale Structure of the Corona and Inner Heliosphere 

Our modeling of the global properties of the solar corona relies on the MHD model to 
describe the interaction of the solar wind with coronal magnetic fields. We initially developed a 
“polytropic model,” in which an adiabatic energy equation with a reduced polytropic index y is 
used (Parker 1963). We have used this model extensively to understand coronal observations. 
Illustrations of the application of this polytropic model to the “Whole Sun Month” (WSM) 
interval, Aug. 10-Sep. 8, 1996, are given by Linker et al. (1999), Breen et al. (1999), Gibson et 
al. (1999), and Posner et al. (1999). 

Although the polytropic model reproduces many large-scale observations of the corona 
with reasonable accuracy (coronal holes, heliospheric-current-sheet crossings, coronagraph 
images), it fails to provide an accurate description of the solar wind. During our present 
program we have improved this aspect of the formulation by modeling in detail the physical 
mechanisms that describe the transport of energy in the corona and solar wind. One- 
dimensional models have demonstrated the importance of including energy transport in 
reproducing spacecraft solar wind measurements (Withbroe 1988; Habbal et al. 1995). 
Accordingly, we have improved the energy equation in our model to include the effects of 
parallel thermal conduction, radiation loss, parameterized coronal heating, and Alfven wave 
acceleration. The model is described in detail by Mikic et al. (1999). The equations we solve 
are described in Appendix A. 

In this improved model we use Spitzer thermal conduction along the magnetic field in the 
low corona, a parameterized collisionless heat flux in the outer corona (Hollweg 1978), a 
coronal heating source, radiation loss, and we solve the equations for the WKB evolution of 
Alfven waves (Jacques 1977), so that we can include the acceleration of the solar wind by high- 
frequency waves. In this way we can reproduce the large-scale properties of the upper 
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chromosphere and transition region, as well as the emission observed from coronal loops. 
While there are many candidate mechanisms, there is no widely accepted theory of corona 
heating (e.g., Parker 1994), so we have chosen to parameterize the coronal heating function. 
Using this approach, we can test which coronal heating mechanisms reproduce observations (by 
comparing simulated X-ray and EUV emission images generated from the model with observed 
images, for example), as discussed in Section 1.1.2. 

In our model we specify the radial magnetic field at the solar surface r = R s (e.g., from 
synoptic magnetic field observations, or from full-disk magnetograms); this field may evolve in 
time (Mikic et al. 1999). We can also match the transverse component of the magnetic field 
when such measurements become available (e.g., from the SOLIS instrument). The boundary 
conditions on the velocity are determined from the characteristic equations along the magnetic 
field. We have found that we can start the model at the top of the chromosphere, at a 
temperature of 20,000 K, allowing us to include the transition region in our calculation (Linker 
etal. 2001). 

Using this formalism we have modeled the structure of the transition region beneath 
magnetic features with different topologies (i.e., open field lines in coronal holes vs. closed field 
lines at neutral lines). We have also self-consistently modeled the formation, support, and 
eruption of prominences, as described below. In the following sections we describe the 
application of this code to several fundamental problems of interest in coronal physics. A 
component of this model that can simulate the inner heliosphere has also been developed, as 
described in Section 1.3.1. 

Even with these improvements to the energy equation, it must be recognized that a single- 
fluid description (inherent in the MHD model) is still a considerable approximation to the state 
of the corona (Habbal et al. 1995; Hansteen & Leer 1995; Hansteen et al. 1997). In particular, 
SOHO observations imply that the electron temperature is considerably lower in the corona 
than the ion temperature. One-dimensional models (e.g., Li, Esser, & Habbal 1997) have 
extended the theory to multiple fluids. Since it is not trivial to include these effects in multi- 
dimensional geometry, we will direct our attention to the single-fluid MHD model first. 
Eventually (not in this proposed effort) we plan to extend our formalism to include multiple 
fluids. 

L L L A Model o f the Solar Wind Including the Tnmsition Region 

In order to study the relationship between coronal heating and the acceleration of the solar 
wind, we used our MHD code to model the solar wind flow in a 2D (axisymmetric) helmet 
streamer configuration (Lionello et al. 2001). We placed our lower radial boundary at the top 
of the chromosphere (at 20,000 K), as described above, to simulate the flow of the solar wind 
through the transition region, into the corona, and into the inner heliosphere. Figure 1 shows 
plasma temperature at increasingly smaller length scales, with magnetic field lines 
superimposed, indicating that the chromosphere extends to higher altitudes in the closed-field 
region above the neutral line, compared to the open-field region. The chromosphere has the 
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Figure 2. The solar wind speed for three different 
sets of parameters: uniform heating, nonuniform 
heating, and nonuniform heating plus Alfven wave 
flux. Nonuniform heating is required to obtain 
slow (equatorial) wind, and Alfven waves are 
needed to accelerate the fast wind. 


Figure 1 . Plasma temperature at increasingly smaller scales with 
superimposed magnetic field lines for the MHD thermodynamic 
model (a,b,c). In panel (d), a Cartesian projection shows the 
structure of the transition region. 
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Figure 3. Comparison of EIT images of an equatorial coronal hole on August 27, 1996 with an MHD model, (a) Open field 

FiT 0 08 S d A aCk fr0m MHD model, (b) EIT 195A image showing the coronal hole, (c) EIT 284A image (d) Simulated 
bl I 284A image using the MHD model with energy transport. 


Predicting Coronal Structure During the 11 August 1999 Eclipse 



Figure 4. Companson between a composite eclipse image created from photographs taken by Fred Espenak in Lake 
hlazar, l urkey with the predicted polarization brightness of the simulated solar corona from our 3D MHD model. The 
projected magnetic field lines from the model are also shown. Terrestrial (geocentric) north is vertically upward The 
eclipse image is copy righted 1999 by Fred Espenak. 1 F 
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lowest altitude near the boundary of the open/closed field region, where the conduction heating 
from the corona is strongest. 

We have compared the results of this “chromospheric” model to those given by a model 
that employs a “radiative energy balance” (REB) boundary condition (as suggested by 
Withbroe 1988). In the REB model, the upper chromosphere and most of the transition region 
are not explicitly included in the calculation (i.e., the boundary temperature is set to, say, 
500,000 K). The plasma density at the boundary is determined by balancing thermal 
conduction from the corona with radiation loss in the transition region, with a small correction 
for enthalpy flow. The advantage of the REB model is that it is not necessary to resolve the 
large gradients in the lower transition region, making for a more efficient numerical simulation 
of the large-scale corona. We have found that the REB model can provide a satisfactory 
approximation (Lionello et al. 2001). However, it is not well suited for investigating 
prominence formation, since it does not contain a chromosphere from which cold and dense 
material can be lifted into the corona. To study prominences we therefore use the 
“chromospheric” model. 

After considerable experimentation, we have produced a model of the solar wind from the 
chromosphere to 1 AU. Figure 2 shows the radial flow speed with three different choices of 
parameters: uniform heating, nonuniform heating, and nonuniform heating with Alfven wave 
pressure. In order to match the observed fast and slow wind velocities, mass fluxes, and 
temperature at 1 AU, the coronal heating scale length needs to be shorter in the streamer belt 
than near the poles, consistent with Withbroe (1988). We have found that an Alfven wave flux 
is needed to accelerate the fast wind to values consistent with observations. Our results agree 
well with generic in situ observations of the fast and slow solar wind, as well as the observed 
properties in the low corona (density and temperature contrasts between the streamer belt and 
coronal holes). 

7.7.2. EUV Emission Ima ges 

The improved energy flow description in our MHD model makes it possible to model 
coronal EUV emission, just as we have previously done for polarization brightness (pB). The 
temperature obtained from the solution can be used to predict the abundance of the coronal iron 
species and produce “simulated” EUV images. The 17lA, 195A, and 284A emission lines 
observed by EIT arise from the excitation of iron ions (Fe IX, Fe XII, and Fe XV). The iron 
population in the corona is especially sensitive to temperature and has been modeled in the 
CHIANTI package (Dere et al. 1997). The emission is dominated by collisional excitation 
from electron impact, and the apparent emission rate R can be written as 

7 i« I n e 2 (e-^ 17 /V7) F(T) dl , 

los 

where the integral is taken along the line of sight. The abundance function for a particular ion 
species, F(T), as well as the energy of the emitted photon, A E, depends on the line being 
observed. 
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To illustrate the idea, we show a comparison we performed for the Whole Sun Month 
period. In Figure 3 we compare simulated emission from our improved 3D MHD model with 
SOHO/EIT images for 27 August 1996. The EIT images show that an equatorial coronal hole 
(the “elephant’s trunk”) corresponds to open magnetic field lines (as predicted by the MHD 
model). Although there is not sufficient resolution in the calculation to reproduce the fine 
structure seen in EIT, the model does show low emission in the vicinity of the coronal hole. We 
can use a similar procedure to develop simulated Yohkoh soft X-ray images. We can also study 
the high-resolution emission on active-region length scales to interpret TRACE EUV emission 
images. This capability is essential to our proposed investigations of active regions and eruptive 
phenomena. 

7-7.3. C omparison with Eclipse Observations 

We have continued our tradition of predicting the structure of the corona prior to total 
solar eclipses. So far we have made four predictions using magnetic field data from the 
previous solar rotation. Our last prediction, of the 1 1 August 1999 total solar eclipse, was our 
most challenging yet, since this eclipse occurred close to solar maximum, when the structure of 
the corona was considerably more complicated than in previous cases. Our prediction was 
posted prior to the eclipse on the World Wide Web (http : / /haven . saic . com). Figure 4 
compares our prediction with an image of the corona taken by Fred Espenak in Turkey, 
suggesting that our model was able to capture the complex magnetic field and streamer structure 
of the solar maximum Sun (Mikic et al. 2000). Our prediction was not as accurate as for 
previous eclipses, which occurred close to solar minimum. We are planning to predict the state 
of the corona during the forthcoming total solar eclipse in June 2001, which will be seen in 
Southern Africa. 

1.2. The Physics of Prominences and Eruptive Phenomena 

7-2.7. Photosnherie Flux Changes and Solar Activity 

Eruptive solar phenomena, such as coronal mass ejections, prominence eruptions, and 
solar flares, are believed to be initiated by the release of energy stored in the coronal magnetic 
field. The mechanism by which this energy is released is not well understood. Motivated by 
observations showing that magnetic flux emergence is associated with filament disappearances 
(Feynman & Martin 1995), we have found that the emergence of new magnetic flux can lead to 
the disruption of coronal magnetic fields and release of magnetic energy. This can occur when 
magnetic flux of opposite polarity emerges in the vicinity of a neutral line, canceling some of the 
existing flux. Recent observations have shown that magnetic flux cancellation is active at 
filament sites (Litvinenko & Martin 1999). Flux cancellation refers to the disappearance of 
magnetic fields of opposite polarity at the neutral line (Martin et al. 1985). We have found that 
flux cancellation near a neutral line, when it does not exceed a threshold value, can lead to the 
formation of stable magnetic flux ropes in sheared 3D arcades (Amari et al. 1999). When the 
amount of flux cancellation exceeds this threshold, the configuration erupts (Amari et al. 2000). 
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Figure 5 illustrates this process. A side view of the magnetic field lines, with the line of sight 
perpendicular to the neutral line, is shown. Initially, the field evolves quasi-statically; at 
t = 80 ta, when the instability threshold is crossed, the structure suddenly erupts upward on an 
ideal MHD time scale. This result is remarkably similar to a simplified analytic model of 
Sturrock et al. (2001), and the triggering mechanism is supported by observations of filament 
eruptions (Wang & Sheeley 1999). Our formalism for incorporating photospheric magnetic 
flux changes into the boundary conditions of MHD simulations is described by Mikic et al. 
(1999) and Linker et al. (2001). 

Flux ropes have long been considered a candidate for explaining prominence support 
(e.g., van Ballegooijen & Martens 1989, 1990; Rust & Kumar 1994; Chen 1996; Aulanier & 
Demoulin 1998). A significant fraction of interplanetary CMEs are “magnetic clouds” (i.e., 
flux ropes, Burlaga 1988; Gosling 1990; Marubashi 1997; Bothmer & Rust 1997). The 
mechanism we have described is similar to that proposed by van Ballegooijen & Martens 
(1989); our calculations show that prominences can form in the lower corona; their eruption can 
disrupt helmet streamers, ejecting CMEs into the solar wind. These results are described in the 
following sections. 

7.2.2. Pwminence ForTnation and Eruption 

Prominences (also called filaments when observed on the solar disk) support cool, dense 
chromospheric material (at ~ ltF K and 10^®— 10^ 1 cm ^) against gravity in the surrounding hot, 
tenuous corona (~ 10^ K and 10^—10^ cnr^). They are observed to lie above magnetic neutral 
lines in the photosphere and near the base of helmet streamers. Prominences have been studied 
for many years, yet the means by which these structures form and are maintained are still not 
understood, nor is their violent eruption. Three main difficulties confront any prospective 
theory attempting to describe the formation and evolution of prominences: (1) Finding a 
magnetic configuration with “dips” (concave upward portions of flux tubes) that can 
gravitationally support the dense material (e.g., van Ballegooijen & Martens 1989, 1990; 
Antiochos et al. 1994; Aulanier & Demoulin 1998; Amari et al. 1996, 1999); (2) understanding 
the mechanism by which chromospheric material is trapped in the dipped field lines and held 
there to form a condensation (Poland & Mariska 1986; Mok et al. 1990; Antiochos & 
Klimchuk 1991; Antiochos et al. 1999a); and (3) elucidating the process that leads to the 
release of magnetic energy and the disruption of these structures (van Ballegooijen & Martens 
1989; Priest & Forbes 1990; Isenberg et al. 1993). Item (3) is closely related to the problem of 
CME initiation (Forbes & Priest 1995; Linker & Mikic 1995; Low 1997; Mikic & Linker 
1997, Wu & Guo 1997; Antiochos et al. 1999b; Lin & Forbes 2000), since these phenomena 
are linked observationally (Hundhausen 1997) and require the release of stored magnetic field 
energy and the opening of previously closed magnetic field regions (e.g., Aly 1984; Sturrock 
1991; Forbes 1992; Mikic & Linker 1994; Antiochos 1998). 

In the theoretical investigations referenced above, the focus was on modeling individual 
aspects of the problem. A complete picture of prominence formation, evolution, and eruption 
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Eruption of an Arcade Caused by Flux Cancellation 


F ig“ re 5. A time sequence showing the dynamic evolution of the magnetic field lines of a highly sheared bipole region 
when new flux of opposite polarity emerges from underneath the photosphere. The new flux begins to emerge at t = 0 At 
t = 80, an instability makes the arcade erupt violently. 
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Figure 6. Evolution of the plasma density (a) and temperature (b) with superimposed magnetic field, in the self-consistent 
thermodynamic MHD model. Flux cancellation underneath the sheared arcade leads to the formation of a cold, dense 
tilament-like structure. It is later ejected when the flux reduction reaches the critical threshold for eruption. In panel (c) we 
show magnetic field lines, colored by temperature (top) and density (bottom). The helical field supports the cold and dense 
prominence material against gravity. 


3D CME Eruption: Magnetic Field Topology 
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Figure 7. Magnetic field topologies arising from a 3D CME computation, (a) Flux rope field lines (viewed from above 
the north pole of the Sun), (b) An isosurface showing the enhanced plasma density associated with the CME. (c) 
V j^ eare<1 ^ nes overla y the erupting structure; these loops are carried out into the solar wind as part of the 
CME. (d) Disconnected field lines are present behind the flux rope. 
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ultimately requires a comprehensive model of all of these processes. This is particularly true if 
the models are required to reproduce observations (e.g., by producing simulated emission that 
can be compared with SOHO or TRACE images). We have used our MHD model with energy 
transport to demonstrate that we can study the problem of prominence formation and eruption 
with a self-consistent approach. We developed a sheared helmet streamer configuration that 
includes the upper chromosphere and transition region (see Section 1.1.1). When flux 
cancellation is imposed in this model, chromospheric material is trapped on helical field lines 
and lifted into a stable configuration in the lower corona, as shown in Figure 6, and described in 
detail by Linker et al. (2001). Observations consistent with this picture have been reported by 
Lites et al. (1995). With further flux cancellation, the entire configuration erupts into the outer 
corona. 

7.2.5. Interaction of Coronal Loops 

High-resolution images from TRACE and SOHO show that interaction of plasma loops 
can lead to dynamical behavior (Pevtsov et al. 1996). The interaction of straight flux tubes, for 
instance, has been studied by Kondrasov et al. (1999), and the energy release was estimated by 
Melrose (1997). We have investigated the dynamic interaction of several two-loop 
configurations which are most likely to release magnetic energy rapidly, for various relative 
orientations between the loops (Mok et al. 2001). We have found that the interaction of loops 
with opposite magnetic helicity tends to be the most violent. The energy released is sufficient to 
explain a small flare (Tang et al. 2000). 

1.2.4. Propagation ojCMEs to_lAU 

In collaboration with D. Odstrcil and V. J. Pizzo of NOAA/SEC, we have studied the 
propagation of a simulated CME (initiated by flux cancellation) to 1 AU. The CME initiation 
and propagation in the inner corona was computed using our (polytropic) coronal MHD model, 
and the subsequent propagation through the heliosphere was carried out using the NOAA/SEC 
model. We have demonstrated that this approach works well (the CME passes smoothly 
between the adjoining boundaries of the calculation without spurious reflection or wave 
generation), and is significantly more efficient than a single calculation from the Sun to 1 AU. 
An interplanetary shock wave forms ahead of the CME, demonstrating that this is a “fast” 
CME. Initial results were presented at the Fall AGU 2000 meeting (Linker et al. 2000); 
additional results were presented at the Spring AGU 2001 meeting (Odstrcil et al. 2001; Riley 
et al. 2001a), and a manuscript is in preparation. 

7.2.5. The Three-Dimensional Structure and Topology, of CM F. a 

Using polytropic MHD computations we have investigated the 3D structure and 
propagation of CMEs initiated by magnetic flux cancellation. As in the 2D case, flux 
cancellation of sheared magnetic fields leads to the formation of a flux rope along the magnetic 
neutral line; in the 3D case the flux rope is anchored at both ends. Figure 7 shows magnetic 
field lines in the flux rope and isosurfaces of the scaled plasma density when the simulated 
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CME reaches the outer corona. These helical field lines remain attached to the Sun; a spacecraft 
intercepting this structure would not only measure field signatures similar to magnetic clouds, 
but would also observe a signature of bidirectional heat flux which often accompanies CMEs 
(Gosling et al. 1987). Figure 7 shows other topological features associated with the eruption, 
including U-shaped field lines behind the magnetic cloud, which may be associated with heat 
flux dropouts (McComas etal. 1989). 

13. Global Modeling of the Inner Heliosphere 

7.5. 1 . An Empirical Heliospheric Model Driven bv the Coronal Magnetic Field Tnpnlnov 

Global models of the inner heliosphere (e.g., Pizzo 1994a,b; Riley et al. 1999; Riley et al. 
2001b,c,d) can provide a necessary contextual basis with which to interpret in situ 
measurements from one or more spacecraft. Other approaches to modeling solar wind 
properties have been described by Usmanov et al. (2000) and Guhathakurta & Sittler (1999). 
As part of our current investigation, we have added a heliospheric component to our 3D MHD 
coronal model, enabling us to model the large-scale structure of the inner heliosphere out to 5 
AU (Riley et al. 2001b,c,d). Since the time step required to advance the solution in the 
heliospheric component of the model is considerably larger than the time step required for the 
coronal solution, we separate the region of space between the solar photosphere and the Earth 
(or beyond) into two parts; the “coronal” region, which includes the region from the 
photosphere up to ~ 3()7? s and the “heliospheric” region, which covers the region between 
30 Rs and 5 AU, where the flow is supersonic and superAlfvenic. We run the two components 
sequentially, with output from the coronal component providing the input into the heliospheric 
component. We have demonstrated that splitting the domain into two separate regions is well 
posed, since we were able to successfully propagate a CME from the coronal model into the 
heliospheric model, as described in Section 1.2.4. As discussed in detail by Riley et al. 
(2001b,c), the heliospheric model currently has an empirical component in the specification of 
the heliospheric boundary conditions. It was developed as an interim solution, so that we can 
immediately begin to model the structure of the solar wind in the inner heliosphere in an 
operational way. In our proposed program we will eliminate this empirical prescription; in the 
future, our heliospheric model will be driven directly using output from the coronal model that 
incorporates a more accurate energy equation. 

To illustrate our results, we summarize the global structure of the heliosphere during the 
WSM interval in Figure 8. The heliospheric current sheet is displayed out to 5 AU, and a 
meridional slice of the radial velocity is shown at an arbitrary longitude. Blue corresponds to 
slowest speeds (~ 350 km/s) and red corresponds to fastest speeds (~ 750 km/s). 
Superimposed is a selection of interplanetary magnetic field lines, as well as the trajectories of 
the WIND and Ulysses spacecraft. The structure portrayed by Figure 8 is consistent with the 
general picture deduced from solar and interplanetary observations during this time period (e.g., 
Riley et al. 1999; Linker et al. 1999; Riley et al. 2001b,c). In Figure 9 we compare our MHD 
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Figure 10. Evolution of several solar parameters during solar 
cycles 21, 22, and the ascending phase of 23, with emphasis 
on the evolution of the HCS during solar cycle 22. The image 
shows the average radial magnetic field from Kitt Peak. 


Figure 11. Comparison of Ulysses solar wind speed 
with the solar wind speed deduced from 3D MHD 
simulations. The model was run for each Carrington 
rotation designated by the vertical bars in the lower 
panel. 
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results directly with Ulysses in situ measurements, showing that the simulation successfully 
reproduces Ulysses’ encounter of an intermediate solar wind stream. 

7. i.2. Solar Cycle Dependence ofjhe Heliosphere 

We have also used this model to explore the evolution of the heliospheric current sheet 
(HCS) during the course of the solar cycle (Riley et al. 2001d). In Figure 10 we present 
several solar parameters, measured over a period of ~ 2-1/2 solar cycles, together with a 
selection of simulation results. The two central panels include data from solar cycles 21, 22, 
and the ascending phase of cycle 23. The lower-central panel shows the average tilt angle 
(black) of the HCS as derived from photospheric measurements using the WSO source-surface 
model, smoothed over 3 Carrington rotations, together with monthly (yearly )-averaged values of 
Sunspot number in red (blue). The upper-central panel is a form of the so-called “butterfly 
diagram” summarizing the longitudinally averaged radial component of the magnetic field, as 
inferred from Kitt Peak magnetograms. Blue indicates inward polarity and red indicates 
outward polarity. 

The HCS for eleven Carrington rotations, spanning solar cycle 22 and covering mid- 1986 
to mid- 1996, is shown above and below the central panels. Inspection of these isosurfaces 
reveals several noteworthy features. First, surrounding solar minimum, the HCS is better 
described as a flat surface with one or more folds in it, in contrast to the sinusoidal picture 
implied from the simple interplanetary extension of a tilted dipole. Second, folds in the HCS 
are typically asymmetric with respect to heliocentric distance: a fold rises more sharply on the 
inner radial side and falls more slowly on the outer radial side. This is a natural consequence of 
the dynamic interaction of the surrounding streams and is particularly effective near solar 
minimum. Riley et al. (200 Id) discuss these results in more detail. 

1-3.3. Understa nding In Situ Solar, Wind. Measurements. 

The heliospheric model has proved to be helpful in interpreting in situ solar wind 
measurements from WIND, ACE, and Ulysses during a number of phases of the solar cycle 
(Riley et al. 2001b,c,d), reproducing many of the large-scale features of the observations. We 
have also used the near-solar-maximum results to speculate on the structure of the high-latitude 
solar wind that Ulysses would encounter during its traversal of the southern and northern solar 
poles in 2000 and 2001, respectively (Riley et al. 2001c). Our results suggested that, due to the 
presence of equatorial coronal holes, the ordered pattern of CIR tilts and their associated 
shocks, which was observed during Ulysses’ initial southward excursion in 1992, would likely 
disappear completely as Ulysses moves toward the South Pole. We also suggested that 
Ulysses would encounter fast streams but would not remain within them for more than a 
fraction of a solar rotation, and that crossings of the HCS would persist up to at least ~ 70° 
heliographic latitude. The most recent measurements from Ulysses (April 2001) are in general 
agreement with these results. 

In Figure 11 we compare Ulysses/SWOOPS wind speed measurements with our 
simulation results, obtained by “flying” the Ulysses trajectory through a sequence of 
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simulations. The solid vertical lines in the lower panel mark the Carrington rotations for which 
we ran simulations. Thus prior to 1999, the comparison is of limited value, since approximately 
1 year is represented by a single Carrington rotation. Nevertheless, this comparison serves to 
illustrate how simulation data can be used to provide a global context for in situ solar wind 
observations. It is apparent that the MHD model reproduces the large-scale changes associated 
with transitions between fast and slow solar wind streams encountered by Ulysses. 
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MHD EQUATIONS USED IN THE GLOBAL CORONAL MODEL 
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MHD EQUATIONS 

For MHD computations to be accurate enough to compare with detailed emission images 
and predict solar wind properties at 1 AU, significant extensions to the usual MHD equations 
are necessary. Following earlier one-dimensional models (Withbroe 1988; Habbal et al. 1995), 
the energy equation in our model to include the effects of parallel thermal conduction, radiation 
loss, parameterized coronal heating, and Alfven wave acceleration, and solves the following 
MHD equations: 



VxbA',, 

(Al) 


VxE.-iS. 

c dt ’ 

(A2) 


E + ^vxB = 77 J , 

(A3) 


^ + V (pv) = 0 , 

(A4) 

^(^ +VVV ) 

= ^JxB-Vp- Vp w + pg + V(vpVv) , 

(A5) 

+ V (pv) = (y - 1 )( - pVv + 5 ) , 

(A 6 ) 

S = - 

V q - n e n p Q(T) + H ch + H d + D , 

(A7) 


where B is the magnetic field, J is the current density, E is the electric field, p, v, p, and Tare the 
plasma mass density, velocity, pressure, and temperature, and the wave pressure p w represents 
the acceleration due to Alfven waves. The gravitational acceleration is g, y = 5/3 is the ratio of 
specific heats, t] is the resistivity, v is the viscosity, H c h is the coronal heating source, D is the 
Alfven wave dissipation term, n e and n p are the electron and proton density, and Q(T) is the 
radiation loss function (Rosner et al. 1978). The term H d = t]fi + vVv:Vv represents heating 
due to viscous and resistive dissipation. In the collisional regime (below ~ 10/? 5 ), the heat flux 
is q = -Kjjbb VT, where b is the unit vector along B, and Kj| = 9 x lO 7 P /2 is the Spitzer value 
of the parallel thermal conductivity. In the collisionless regime (beyond ~ 10 R s ), the heat flux is 
given by q = cxn e kT\, where a is a parameter (Hollweg 1978). Since it is presently not known 
in detail what heats the solar corona, the coronal heating source H c \^ is a parameterized function. 
A simple form is 

H ch = H 0 (6) exp [- (r - R 5 )/A( 0 )] , (A 8 ) 

where H 0 (£f) expresses the latitudinal variation of the volumetric heating, and A(0) expresses the 
latitudinal variation of the heating function scale length. [In practice, this function can be 
tailored to match the scaling of different coronal heating models). Note that the simplified 
polytropic model is obtained by setting S = 0 in Eq. (A 6 ), p w = 0 in Eq. (A5), and y = 1.05. 


18 



f 


SAIC-01/801 0:APPAT-277 
June 4, 2001 

Since the acceleration of the solar wind by Alfven waves occurs on a spatial and time scale 
that is below the spatial and time resolution of our global numerical model, the wave pressure 
Pw is evolved using an equation for the time-space averaged Alfven wave energy density e 
(Jacques 1977), 

ds 

37 + V F = v-Vp w -D , (A9) 

where F = (- v + \ A )s is the Alfven wave energy flux, v A = BI^Anp is the Alfven speed, and 
Pw = 2 e - The Alfven wave velocity is V 4 = ±bv^; in a multi-dimensional implementation, it is 
necessary to transport two Alfven wave fields (waves parallel and antiparallel to B), which are 
combined to give s. The Alfven wave energy density e is related to the space-time average of the 
fluctuating component of the magnetic field SB by £=<SB 2 >/4jc. The dissipation term D 
expresses the nonlinear dissipation of Alfven waves in interplanetary space and is modeled 
phenomenologically (Hollweg 1978). 

We have developed three-dimensional codes to solve MHD equations (A1)-(A9) in 
spherical coordinates (r ,G,(pi) and equations (A1)-(A7) (without Alfven wave acceleration) in 
Cartesian coordinates. The spherical code have been used extensively to model CMEs (Mikic 
& Linker 1994, 1997; Linker & Mikic 1995, 1997; Linker et al. 2001), coronal structure 
(Miki 6 & Linker 1996; Linker et al. 1999; Mikic et al. 1999, 2000; Lionello et al. 2001), and 
heliospheric structure (Riley et al. 2001b,c). The Cartesian code uses a very similar algorithm 
and has been used to study active regions (Mikic & McClymont 1994; Mikic et al. 1996; Jiao 
et al. 1997; McClymont et al. 1997; Mok et al. 2001). 
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